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Abstract 

The simulation of runoff from a 
Himalayan Glacier basin using 
an Artificial Neural Network 
(ANN) is presented. The 
performance of the ANN model is 
found to be superior to the 
Energy Balance Model and the 
Multiple Regression model. The 
RMS Error is used as the figure 
of merit for judging the 
performance of the three 
models, and the RMS Error for 
the ANN model is the least of 
the three models. The ANN is 
faster in learning and exhibits 
excellent system generalization 
characteristics . 

Introduction 

Central Asia, the main regime 
of mountain glaciation in the 
Himalayas, is notable for its 
diversity of hydrological and 
meteorological conditions. As 
the largest ice-sheet outside 
the polar regions, it nourishes 
some of the largest river 
systems, like the Indus, the 
Ganges, etc. These mountain 


rivers have as their source 
the liquid run-off from melting 
snow and ice as well as liquid 
precipitation into their 
drainage basins. However, for 
most of the river systems in 
the Central Asia, the water 
from melting is the main source 
of nourishment. 

The study of the glaciated 
mountainous areas presents many 
difficulties, due to both lack 
of information on the water 
balance elements and the 
great variety and complexity of 
natural conditions of these 
areas. The structure of natural 
zones of ablation ( and 
consequently the water balance) 
is determined by geographical 
position and the peculiarities 
of orography, the absolute and 
relative height of mountains, 
the circulation processes and 
radiation regime, the 
exposition of slopes, the 
geological structure, etc. [7] 

There exist many empirical 
methods of simulating and 
predicting glacier melt runoff, 
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including time series modeling 
and actual physical parametric 
modeling. One of the most 
widely studied and viable 
method is the Energy-Balance 
Model. This method has been 

advocated by various 

researchers for operational 

runoff forecasting from 

glacierized basins. [3]/ [4]/ 

[ 10 ] 

The Study Area 

The Chhota Shigri glacier is a 
medium size valley glacier 
situated in the Lahaul and 

Spitti District of Himachal 
Pradesh (India) . The glacier 
lies between latitudes 77-29N 
and 77-33N and longitudes 32- 
11E to 32-18E, covering an 
area of about 19.39 sq. km. 
Geologically this area is 
confined to the Rohtang 
gneissic complex of the Central 
Crystalline of the higher 
Himalaya. Fig. 1 (a) shows 

various features of the glacier 
and Figure 1 (b) shows relevant 
hydrogeomorphological 
parameters. [1] 

The Energy Balance Model 

Marks and Dozier (1992) have 
carefully analyzed the various 
types of heat fluxes and their 
relative influence over the 
snow melt run-off process in a 
small alpine basin. They have 
clearly demonstrated the 

dominance of solar radiation 
heat flux and two other 
important fluxes, viz. the 
sensible heat flux, and the 
latent heat flux in the runoff 
calculations. This logic can be 
applied to other glaciated 
areas as well. Since the solar 
radiation is easily measured 
and the heat flux is easily 



Figure 1 (a). Glacier Location Map 


Total Area of the Glacier 

19.39 sq. km. 

Accumulation Area 

18.15 sq. km. 

Average Height of Acc. Area 

5010 Mts. 

Ablation Area 

1.24 sq. km. 

Average Height of Abl. Area 

4133 Mts. 

Accumulation Area Ratio 

0.91 

Maximum Glacier Height 

5595 Mts. 

Snout Height 

3840 Mts. 

Orientation 

NNE 

Average Gradient 

176 Mts. /km. 

Length of Glacier 

10.00 km. 

(adopted from Dhanju and Buch, 1989) 


Figure 1 (b). Hydrogeomorphological Parameters 
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Run-off in Cusec. 



Observed Runoff 

Figure 2. Hydrograph of Chhota Shigri Stream 


modeled for a given topography, 
it should be possible to 
reliably model the runoff from 
glaciated areas. 

The energy budget for a melting 
glacier can be expressed as 

Qm = Qn + Qh + Qe + Qr ( 1 ) 


observed during the period of 
this study.) 

The net radiation flux is 
calculated as 

Qn = Qsw + Qlw (2) 

Where 


Where 

Qm = Energy flux actually 
available for melting the ice. 

Qn = Net radiation flux ( in 
Watts / Sq. Meter ) 

Qh = Sensible heat flux ( in 
Watts Sq. Meter ) 

Qe = Latent heat flux ( in 
Watts / Sq. Meter ) 


Qsw = Shortwave radiation 
balance measured at the 
observatory, by net radiometer. 

Qlw = Long wave radiation 
budget at the surface, modeled 
using the algorithm suggested 
by [13]. 

The turbulent fluxes, viz. The 
Qh and Qe were determined using 
the standard aerodynamic 
formulae and the roughness 
length was assumed to be 
constant at 5mm for this 
experiment. 


Qr = Sensible heat flux from 
rain. ( This term is negligible 
in our case as no rain was 
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Thr** Input Unit* 

Figure 3. Neural Network Architecture 
used lor Glacier Run-oil Simulation 


The hydrometeorological data 
and the net radiation flux 
data, along with corresponding 
runoff, were collected in situ 
during the Chhota Shigri 
Glacier Expedition in July - 
August, 1988. The approach to 
this glacier is possible only 
during the ablation season, as 
the high mountain passes remain 
snow-covered for the rest of 
the period. This is the major 
bottle-neck in collecting 
long-term melt-discharge 

records for these areas. The 
Hydrograph of the Chhota 

Shigri melt-stream is shown in 
Figure 2 . 

The data are sparse and call 
for a technique that can 
generalize the latent system 
configuration to improve their 
applicability to runoff 

calculations. One such recent 
technique that possesses 

excellent generalization 

characteristics in addition to 
very good learning capability 
is the Artificial Neural 
Network technique. 

Modeling Physical Systems and 
Processes using Artificial 
Neural Networks 


Many of our research efforts in 
the field of Remote Sensing 
remain directed towards 

understanding and developing 
mathematical models that can 
simulate dynamic physical 
processes in oceans, 

atomosphere and the environment 
in general. The spacecrft-based 
sensors provide much-needed 
details about various 

parameters that govern these 
processes. However, due to 
limited knowledge of the 

physical processes in the 
environment, and the inherent 
noise in many geophysical data, 
environmental systems often 
cannot be accurately 

represented through numeric 
values describing their 

physical properties and 

interactions, but rather are 
subjected to categorisation 
into broad classes. [9] 


The modelers have been 

investigating many different 
techniques to exactly simulate 
the processes from the previous 
knowledge, available in terms 
of historic data (very good 
survey of recent tachniques is 
provided by Tsonis and Eisner, 
1989) . Many efforts to model 
dynamic systems have been 
undermined by the system's 
inherent chaoticity in the 
data. One solution is to model 
the temporal variability of a 
parameter rather than the 
parameter per se. [2] 

Neural Network Architecture 

used for Runoff Simulation 


Recently, there has been a 
great surge in the application 
of multi-layer-feed-forward 
Artificial Neural Networks to 
diverse system identification 
problems [2], [11]/ [5]. In 
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Figure 4. Comparative Model Performance 


training is complete. 

The standard method for 
measuring the generalisation 
characteristics of a given 
model is called the method of 
'Cross Validation'. The method 
splits the data set into two 
subsets, viz. the training data 
set and the testing data set. 
The learning is performed using 
the training data set and the 
network performance is 
evaluated using the test data 
set. Unless the training set 
is large enough, the 
performance of the network on 
training data is not likely to 
be an accurate measure of its 
performance on future or 
unknown data. 

To achieve a statistically 
significant result, several 
independent splits are required 
and the average of the results 
is accepted as the overall 
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this study, the 

Fast-Back-Error-Propagation 
algorithm has been used. [6] 

The network for run-off 

prediction comprises 3 layers, 
viz. the input layer, the 
hidden layer and the output 
layer. The input layer consists 
three input neurons 

corresponding to three heat 
fluxes viz. the Solar Radiation 
Heat Flux, the Sensible Heat 
Flux and the Latent Heat Flux. 
The hidden layer has 30 

nodes, and the output layer 
consists of one neuron 
corresponding to the modeled 
run-off value. The schematic 
diagram -is shown in Figure 3 . 

Generalisation Performance of 
ANN Model 

Generalisation is a measure of 
how well the network performs 
on the actual problem once the 





training data set. [5] 


PSE = MSE+ (2*Nw/Np) * (r2) (3) 


Model RMS Error 


Energy Balance Model 
Multiple Regression Model 
ANN Model ( 1391 Iterations ) 


Figure 5. Table Showing Modeling Errors 


performance of the network. 
While the cross validation 
method is a widely accepted 
method, it is extremely time- 
consuming in the case of ANN as 
lenghty training times are 
required for each independent 
partition of the data set. 


Where 

PSE = Pridicted Squared Error. 

( For Future ) 

MSE = Mean Squared Error of the 
training data set. 

Nw = Number of free parametres 
in the model (120 in our 
case ) . 

Np = Number of training 

patterns for the model (24 in 
our case ) . 

(r2) = Variance of the noise 

determined by the formula : 


In the case of the Himalayan 
glacier basins, 

hydrometeorological 
collected in situ are 
as the mountain passes 
to the glacier snout 
open for a very 

period towards the 
summer. Hence, 
sets have to 
meaningful training. Thus there 
is a need for alternate 
reliable methods for 

predicting the generalisation 
performance of the network 
without having a test data set. 


the 
data 
sparce , 
leading 
remain 
short 
end of 
all of the data 
be used for 


(r2) = abs[ (p/ (p-Nw) ) *MSE] (4) 
Where 

p = Number of patterns used for 
training the network. 

The MSE in the ANN model for 
the Chhota Shigri Runoff Model 
is estimated at 0.0036. Hence 
by substituting the values in 
eq. 4 and eq. 3 we get the PSE 
as 0.0126, which is quite 
acceptable as model error 
considering the paucity of 
data . 


An alternative technique which 
requires far fewer computations 
is called the 'Predicted 
Squared Error' technique. This 
technique relies on statistical 
methods to derive an expression 
for the generalisation 
performance of a system as a 
function of its performance on 
the training data set, the 
number of free parameters in 
the system and the size of the 


Conclusions 

Figure 4 depicts the resultant 
calculated runoff from the 
energy balance model against 
two other models, viz. the 
Multiple-Regression Model, the 
Energy Balance Model and the 
ANN model. (See Figure 5) 

It can be clearly seen from the 
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comparative performance of the 
three different models that the 
trained ANN model assumes the 
flow values that are closest to 
in situ measured values. 
(The Correlation Coefficient 
between the measured runoff and 
the runoff simulated by ANN 
model is of the order of 
0.9998) Currently, we are 
investigating the application 
of different ANN models to 
diverse hydrological and 
meteorological simulation and 
forecasting problems. 

The technique of system 
modeling with ANN holds very 
good promise. However, it 
requires rigorous research and 
simulations before any 
operational model can be 
developed . 
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